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Abstract 

We consider symmetry properties of an integro-difTerential multidimensional Gross-Pitaevskii 
equation with nonlocal cubic nonlinearity in the context of symmetry analysis using the formalism of 
semiclassical asymptotics. This yields a semiclassically reduced nonlocal Gross-Pitaevskii equation 
which determines the principal term of a semiclassical asymptotic solution and can be referred to 
as a nearly linear equation. 

Our main result is an approach which allows one to construct a class of symmetry operators 
for the reduced Gross-Pitaevskii equation. These symmetry operators are determined by linear 
relations including intertwining operators and additional algebraic conditions. 

The basic ideas are illustrated by the examples of a ID reduced Gross-Pitaevskii equation. 
The symmetry operators are found explicitly, and the corresponding families of exact solutions 
generated by the symmetry operators are obtained. 

PACS: 02.30.Jr Partial differential equations; 03.65.Sq Semiclassical theories and applications; 
ll.lO.Lm Nonlinear or nonlocal theories and models; 11.30.-j Symmetry and conservation laws 

MSG: 35Q Equations of mathematical physics and other areas of application; 35Q55 NLS-like 
equations (nonlinear Schrodinger) ; 76M60 Symmetry analysis, Lie group and algebra methods; 
81Q20 Semiclassical techniques, including WKB and Maslov methods 



1 Introduction 

Symmetry operators which, by definition, leave the set of solutions of an equation invariant, 
are of essential importance in symmetry analysis of nonlinear partial differential equations (PDEs) . 
The obvious use of symmetry operators of the equation is generation of new solutions from a known 
solution. 

However, no general approaches are known for regular calculation of symmetry operators for 
nonlinear equations due to the equations that determine symmetry operators are nonlinear operator 
equations. Solving them is a complicated mathematical problem which requires special techniques 
not yet developed. 

In addition, to solve the determining equations, we have to define the structure of symmetry 
operators consistent with the determining equations, but there are no recipes for choosing such a 
structure. As a consequence, finding the symmetry operators for nonlinear equations in general 
problem statement is an unrealistic task. Therefore, the Lie group properties of PDEs and related 
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structures are studied in symmetry analysis of differential equations [IHT] where the main tool is 
the symmetry that is generator of the Lie group of symmetry operators. 

Modification and further development of the classical Lie group techniques for nonlocal (integro- 
differential) equations has a long history (see e.g. [8l-[12]). The idea of current research in this area 
can be formed on the works [13H2T]. 

The symmetries of PDEs are effectively calculated from linear determining equations and widely 
used in finding of classes of particular solutions, conservation laws, etc. 

At the same time, the construction of a Lie group on a given symmetry by integrating of the 
Lie system can be performed explicitly only for a special kind of symmetries. Having an algorithm 
of finding of families of symmetry operators we obtain a possibility to construct Lie groups of 
symmetry operators and consider structure of the correspondent symmetries. 

Note that for linear PDEs, symmetry operators are effectively found from linear determining 
equations and are widely used in applications such as quantum mechanics (see, e.g., pr71l22|j23] and 
references therein). This fact leads to the idea to look for a special class of nonlinear equations for 
which symmetry operators could be calculated using the methods applicable to linear equations. 

An example of such a class of nonlinear integro-differential equations (IDEs) with partial deriva- 
tives is considered here. We call equations of this class the nearly linear equations. The symmetry 
operators for them can be found by solving linear operator equations (similar to linear PDEs) and 
additional algebraic equations. 

We consider a generalized multidimensional integro-differential Gross-Pitaevskii equation (GPE) 
with partial derivatives and a nonlocal cubic-nonlinear interaction term of general form. The 
WKB-Maslov method of semiclassical asymptotics PHHS] is used to obtain a reduced GPE from 
the original GPE. The reduced GPE is quadratic in spatial coordinates and derivatives and con- 
tains a nonlocal cubic-nonlinear interaction term of special form. This equation belongs to the 
class of nearly linear equations and determines the principal term of the semiclassical asymptotic 
solution. 

The main result of our work is an approach developed for finding of symmetry operators for 
the reduced GPE by solving linear operator equations. 

The general way for using this approach is illustrated by the example of a one-dimensional 
reduced GPE whose symmetry operators are found explicitly. With the use of symmetry operators 
obtained, two families of exact solutions are generated for the reduced GPE. 

In the following section, an integro-differential Gross-Pitaevskii equation is introduced and its 
semiclassical reduction is presented. A method for integrating the reduced GPE is described. We 
obtain the foundation stones of the method, namely a consistent system and a linear equation 
associated with the reduced GPE. 

In section 2 we propose an approach for finding the class of symmetry operators of the reduced 
GPE by constructing intertwining operators. 

In section 3 the general ideas are illustrated by the example of a one-dimensional GPE of 
special type. The symmetry operators for this equation are found explicitly and two families of 
exact solutions are generated with the use of the operators obtained. 

2 The nonlocal Gross— Pitaevskii equation and the consis- 
tent system 

The Gross-Pitaevskii equation and its modifications are widely used in studying coherent mat- 
ter waves in the Bose-Einstein condensates (BECs) . Recent extensions to the BEC studies 
involve long-range effects in the condensates described by a generalized GPE containing integral 
terms responsible for nonlocal interactions. We refer to equations of this class as nonlocal GPE 
(which are also known as the Hartree-type equations in the mathematical literature). The nonlo- 
cal BEC models, when used in the theory of BECs with attractive particle interactions, may keep 
the condensate wave function from collapse and stabilize the solutions in higher dimensions (see, 
e.g., [5DH32| . a review ^33,, and references therein). Nonlocal GPEs also serve as basic equations 
in models describing many-particle quantum systems, nonlinear optics phenomena [34], collective 
soliton excitations in atomic chains [3F, etc. 





Let us write down the nonlocal Gross-Pitaevskii equation as 



F(*)(f,t) = {~ihdt + Hit) + xV{^){t)}-^{x,t) = 0, 
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(1) 



V{^){t) = V{^){z,t)^ J dy^*{y,t)V{z,w,t)^{y,t), 



(2) 



where ^{x, t) is a smooth complex scalar function belonging to a complex Schwartz space S in the 
space variable a; € R" at each point in time t. In our consideration, we solve the Cauchy problem 

^ix,t)l^^=i:{x) (3) 

with the initial (at t = s) function tp{x) € §. 

The space § is taken to ensure existence of the Euclidean norm = \/JW(J)^W(Jy) , of the 

moments of i), and of the integral in ©. Here, = J dx^*{x,t)'ii{x,t) denotes 

R" 

the Hermitian inner product of functions ^{x,t) and 'i/{x,t) of the space S; ^*{x,t) is complex 
conjugate to ^{x,t). 

The linear operators H{t) = H{z,t) and V{z,w,t) in ([ij are the Hermitian Weyl-ordered 
functions |36j of time t and of the noncommuting operators 

z = {p,x) ~ {~ihd/dx,x), (4) 
w = {-ihd/dy, y), x,y€R", 

with commutators 



[zk,Zj]^ = [wk,Wj]- ^ ihJkj, [zk,Wj]^^0, fc,j = l,2n, (5) 
where [A, B]^ = AB - BA, J = \\Jkj\\ 2nx2n is a Unit symplectic matrix: J = ( rr ^] , and 



2nx2n 



I 

I = Inxn is an ?T. X n identity matrix. We also use shortcut symbols for derivatives: dt ~ d/dt, 
dg = d/dx. 

From Eq. ([T|) it follows immediately that the squared norm of a solution '^{x,t) of Ec^. ^ is 
invariable with time, = = HV'IP- 

In the multidimensional case, the GPE ([iJ with variable coefHcients of general form is non- 
integrable by known methods, such as the Inverse Scattering Transform [37] ■ Therefore, analytical 
solutions to this equation can be constructed only approximately. An effective approach to con- 
structing such asymptotic solutions is the method of semiclassical asymptotics at ?i — > 0. Thus, 
the theory of canonical operator with real phases was developed by Maslov [SHllSn] for solving of 
the Cauchy problem for self-consistent field equations. Spectral problems with singular potentials 
were considered in [lUlSI] (see also [12H45j ) and soliton-like solutions of the Hartree-type equa- 
tion with potentials of special form were constructed [46] . A specific and attractive feature of the 
nonlocal GPE ([1]) is that, in the semiclassical approximation, the input GPE is reduced to an 
equation containing nonlocal terms which can be expressed as a finite number of moments of the 
unknown function '^{x,t). The reduced equation can be considered as a nearly linear equations. 
The concept of a nearly linear equation implies that among the solutions of a nonlinear equation 
there is a subset of solutions regularly depending on a nonlinearity parameter 1361 . 

According to the WKB-Maslov complex germ method [24l[25] applied to a class of trajectory- 
concentrated functions, the construction of semiclassical asymptotics for Eq. ([1} is reduced to an 
auxiliary problem for an associated linear Schrodinger equations. 

A solution of the Cauchy problem ([!]) and ([3]), asymptotic in a formal small parameter h 
{h 0), was constructed accurate to 0{h^/^) where N is any natural number jJT]. The leading 
term of the asymptotics is found by reducing the GPE ([iJ to a GPE with a quadratic nonlocal 
operator 

{-iHdt + Hg{z,t) + K J dy^*{y,t)V,{z,w,t)<I'iy,t)}^ix,t)=0 (6) 

R" 

where the linear operators Hq{z,t) and Vq{z,w^t) are Hermitian and quadratic in z, w: 

Hq{z,t) = ]^{z,'H.S)z) + {H.{t),z), (7) 
Vq{z,w,t) = ^{z,W^4t)z) + {z,W:,^it)w) + ^{w,W^n^it)w). (8) 



Here, V-zzit), Wzz{t)^ Wzw{t), and W.uiw{i) are 2n x 2n matrices; V-zit) is a 2n vector; the angle 

n 

brackets (., .) denote the Euchdean inner product of vectors: (p, x) — ^ Pj^ji Pt^ £ M"; {z, w) = 

j=i 

Yl'j=i ZjWji z,w Q M^". We shall call equation ([5]) a reduced Gross-Pitaevskii equation. 

The reduced GPE given by Eqs. (HJ, (O, ([8]) can be integrated explicitly [491,150' and it possess 
fairly rich symmetries. Analysis of these symmetries can provide a wealth of information about 
the equation and its solutions. 

As Eq. ([6|) contains a nonlocal nonlinear term, its symmetry properties are of special interest 
in the symmetry analysis of partial derivatives equations. The matter is that the application of 
the standard methods of symmetry analysis [IH31[51[n], developed basically for PDEs, to equations 
different in structure from PDEs presents a number of difficulties: For instance, there are no 
regular rules for choosing an appropriate structure of symmetries for non-differential equations. 
The reduced equation ([6]) allows one to avoid this problem, as the symmetry properties of this 
equation are closely related to the symmetry of the linear equation associated with the input 
nonlinear equation. 

The key factor in analyzing the symmetries of the nonlinear equation F('^){x,t) = is the 
symmetry operator A that makes the set of solutions of the equation invariant (see, e.g., [7ll22j): 

= ^ = 0. (9) 

Generally, it is impossible to find effectively a symmetry operator A for a given nonlinear operator 
F by solving the nonlinear operator equation ([9|) . This situation is resolved in the group analysis of 
differential equations [IHl] where a symmetry a (generator of a Lie group of symmetry operators) 
is the main object of analysis. 

The symmetries are determined by the linear operator equation 

F(*)(f,t) = ^ /"((T*)(.f,i) = 0. (10) 

Here, F'(^') is the Freshet derivative of F calculated for ^ff. For the linear operator F, we have 
F' = F and the symmetry operators being the same as the symmetries. 

We assign the GPE ([6]) to the class of nearly linear equations following the definition given 
elsewhere [48]: A nearly linear equation determining a function ^ has the form of a linear partial 
differential equation with coefficients depending on the moments of the function ^. This type of 
equation can be associated with a consistent system which includes a system of ordinary differential 
equations (ODEs) describing the evolution of the moments and the reduced GPE. 

Using the GPE as an example, we can see that the class of symmetry operators for nearly linear 
equations can be found by solving the respective determining linear operator equations. In this 
sense, the symmetry properties of nearly linear equations are similar in many respects to those of 
linear equations. 

Let us consider briefiy a method for solving the Cauchy problem Q for the reduced GPE ^ 
following the scheme descibed elsewhere [37] • 

We denote the Weyl-ordered symbol of an operator A{t) = A{z,t) by A{z,t) and define the 
expectation value for A{t) over the state ^'(a?, <) as 

= = ^/ dx^*{x,t)Am{x,t). 

As II^'IP does not depend on time, we have from ([6]),(l7]), ([8]): 

= ^ / dx^*ix,t){^ + ^[H,iz,t)J{t}].+ 

R" 

+ dy^*{y,t)[Vg{z,w,t),A{t)]^^ix,t)Y (11) 

R" 

where A^s,(t) = dA^{t)/dt. 

We call Eq. similar to the linear Schrodinger equation in quantum mechanics (>f = in 

Eq. H])), the Ehrtenfest equation for GPE 
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Let z^i,{t) — (z,p;(t)) and A^^\t) ~ (A|^].;(t)) denote the expectation values over for 
the operators 



k,l = 1, 2n, 



(12) 



(2) 

respectively. Here, Azi — zi — We call the first moments and AJj, (t) do the second 

centered moments of i). 

From ([6]), ([7]), ([8]), and (jlll) we immediately obtain a dynamical system in matrix notation: 



'i* = J{H,(i) + + ?c{W,,{t) + W,^,{t))]z^}, 



(13) 



The functions g = 0*(i) describe phase orbits in the phase space of system (|T3)) . 
Then the Cauchy problem ^ for the reduced GPE © can be written equivalently as 



We call Eqs. ([T^ i/ie Hamilton-Ehrenfest system (HES) of the second order for GPE ©. The 
HES is of the second order, as Eqs. (|13p contain the first and second moments. 

For brevity, we use a shorthand notation for the total set of the first and second moments of 
*(f,0: 

Q^{t) = {z^{t),A^^\t)). (14) 



(15) 

(16) 
(17) 
(18) 



'-.{t,Q^{t))-^{x,t) = { 



ihdt + H, 



,(t,fl*(t))}*(f,t) 



0, 



Hq{t,Q^{t)) = ^{z,n,,{t)z) + {n,{t),z) + ^{z,W,,{t)z) + 

+ ^{z^,{t),W^,^,{t)z^{t)) + }c{z,W,^{t)z^,{t)) + ^Sp 

fl*(t) = r(t,0*(i)). 



,{t) 



Equation pT)) is a concise form of HES (IT3)) . and r(t,g,i,(t)) designates the rhs of ^3]) . 

We call GPE ([T5|) and the correspondent HES (fTT]) the consistent system for GPE 

The GPE (fT5|) can be assigned to the class of nearly linear equations [15] , as the operator 
of the GPE (fT5|) is a linear partial differential operator with coefficients depending only on the first 
and second moments g*(t). 

The consistent system (|17l) allows us to reduce the Cauchy problem for the GPE (fT5|) to 
the Cauchy problem for a linear PDE owing to which the Cauchy problem for HES (|17p can 
be solved independently of Eq. (|15p . 

Let 

fl(i,C)= (z(i,C),A(2)(t,C)) (19) 

be the general solution of HES p7p and C ~ (Ci,C2,..., Cat) denote the set of integration 
constants. 

Consider a linear PDE with coefficients depending on parameters C: 



L{t,C)^{x,t,C) = { - ihdt + Hq{t,C)}<S>{x,t,C) = 0, 



(20) 



where 



Hg{t,C) = ^{z,HUt)z) + {n,{t),z) + ^{z,W,,{t)z) + 
+k{z, WUt)Z{t, C)) + ^{Zit, C),W^Ut)Z{t, C)) + I Sp 



iy„^(0A(2)(i,C) . (21) 



The operator Hq{t, C) (l?T|) of Eq. (I^Ul) is obtained from p^ where the general solution g(t, C) of 
the HES PT)) stands for the moments g*(t). 

We call Eq. (^0)) t/ie Zmear associated equation (LAE) for GPE p^ . 
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Let ^{x,t,C[tp]) denote the solution of the Cauchy problem for LAE ^U\\ with the initial 
condition 

<i>(x,i,C[V']) (22) 

where integration constants C are replaced by the functionals C = C[-0] determined from the 
algebraic conditions 

fl(i,C) (23) 

t—s 

Then the solution of the Cauchy problem US]), (HH) for the GPE (see [TfllSD] for details) is 

^ix,t)^^x,t,C[ij]). (24) 

Define C[^'](t) by the condition 

0(f,C[*](i))=0*(i). (25) 
From the uniqueness of the solution of the Cauchy problem for the HES (|18p it follows that [50] 

Qit,C[m))^Bit,C[^]). (26) 

and, hence, 

i.e., the functionals C[5'](i) are the integrals of Eq. 
Also, we have 

fl(i,C[7^])=0^(t), (28) 

where 0^(i) is the solution of HES P7|) with the initial condition (P^ . 

From ([26]) we see, in particular, that the parameters C[5'](i) related to 'i'{x,t) do not depend 
on time, i.e. C[*](i) = C[i/']. 

When the symbols H{z, t) and V{z^ w, t) of the operators ((!]) and ([2]) are not quadratic functions 
of z and w, the operator H{t, C) in ([2T|) can be presented as 



m, c) . nz, t) + ^± ^ ""lr''^ A(H)(t, c) 



(29) 



z{t,C) 



Here, £ Z^" is a multi-index 



1/ = (I'l, . . . ,Z/2«), l^k>0, fc=l,2n, 1^1 + . . . + l^2n, 

z/! = (z/i!)---(i/2„!), a,..,. = TTir, TTj;^, 



dw^ dwi^ . . . dw. 



2n 



and A[j'^^^ is a |z^|-order centered moment of the solution ^'(af, t). In other words. Ay"'-* is the 
expectation value of the operator A';j'^^\t) = A^i!'^^\z,t) with the Weyl symbol A^'^^\z,t) = (z — 
z^,{t)Y. The associated linear equation (pO)) is in general not exactly integrable, and to find its 
approximate solutions, we use the semiclassical asymptotic formalism [47(l52j. Analysis of the 
LAE of general form involves a great number of additional technical issues associated with the 
semiclassical approximation that requires a separate study. That is why to illustrate the main 
ideas of the proposed approach, we should limit our discussion to the case of a quadratic operator 
for which Eq. pO| is integrable. 

Let us now turn to the construction of symmetry operators for Eq. ([6]). To do this, we use an 
operator intertwining a pair of LAEs of the form ([20]) . although the symmetry operators can also 
be found by another methods (see also (51j ). 
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3 Intertwining operator and symmetry operators 

According to the definition (j9|), the nonhnear symmetry operator A{t) maps any solution ^(x, t) 
of Eq. (|15p into its another solution: 



= (30) 
For a — A{t)\^_^ and tpix) given by (O, we can set 

V'a(x)=a^(x) = 4'AW|,^,, (31) 

and g^^ are the first and second moments of ipaix) similar to ([T^ . 

From the solution of the Cauchy problem for the HES (|T8| with the initial condition g*(t) = 
g^^, by analogy with (|28|) . we have 

Q{t,C[^ba])=Q4,At)- (32) 
According to ([M]), the solutions *(x, t) and '^'A{x,t) of the GPE ([11]) are found as 

*(f,t) = $(f,t,C)|^^^j^j (33) 

and 

vI'^(f,i) = $(x,t,C')|c,^c'[^„]' (34) 

where t, C) and $(a?, C) are the solutions of two LAEs of the form pop with two different 
sets of integration constants C and C, respectively. 

To construct the symmetry operator A{t), we relate the functions ^{x, t, C'[ipa]) and $(a;, t, C[ip]) 
by a linear operator M{t, s, C , C) which intertwines the two operators L{t, C) and L{t, C) in (|20l) : 

C')M{t, s, C, C) = i?(t, s, C, C)L(t, C). (35) 

Here, the linear operator i?(i, s, C',C) is a Lagrangian multiplier, and the initial condition is 
M{t,s,a,C)\t=s=a. 

From (|35|) we have that $(a;, C) = M{t, s, C, C)<I'(af, i, C) for two arbitrary sets of constants 
C and C, and this is especially true for ^{x, t, C'[ipa]) and $(a?, t, C[ip]) with the constants C'[V'a] 
and C[V']. 

To find the operator M {t, s, C, C), we consider a linear intertwining operator s, C, C) for 
L(t, C) and L(t, C) satisfying the conditions 

L{t, C')T){t, s, C, C) = Vit, s, C, C)L{t, C), (36) 
P(i,s,C',C) =1 (37) 

We call I?(t, s, C, C) t/ie fundamental intertwining operator for L(t, C) and C). With the 
use of T>{t, s, C, C), the operator M{t, s, C, C) of (|35|) can be presented as 

if(i, s, C, C) = 25(t, s, C, C)i3(t, C). (38) 

Here, B{t, C) is the linear symmetry operator of LAE ([20|) satisfying the conditions 

[L(t,C),B(t,C)] =0, B{t,C)\t=s=a. (39) 

Hence, given the operator T>{t, s, C, C) of (|36|) and the family S of linear symmetry operators 
of LAE (p9| we can construct the family of nonlinear symmetry operators for the GPE ([T]): 

(i(t)*)(x,t) =25(t,s,C'[a^A],C[*](t))B(t,C[*](t))*(x,0- (40) 

Here, C'[aV'] and C[vp](= C[V']) arc found from and ([25]), respectively, and B e B. 

Note that the symmetry operator A{t) from Eq. (|40| is nonlinear since the operators M and 
i? depend on the parameters C which are functionals of the function ^1/. 
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To find the fundamental intertwining operator s, C, C), we introduce a function 4'{x, t, C) 
by the conditions 

$(f , t, C) = k{x, t, s, C)(j){x, t, C), (41) 
k{x, t, s, C) = exp[-(X(t, C), V)] cxp [^[S{t, C) + (P(i, C), f)]}, (42) 

where $(f , C) is a solution of Eq. (gDl), the vector z = C) = (P(t, C),X{t, C)) satisfies Eq. 
p7)) . and 5(i, C) is a smooth function to be determined. 
For i, C) we have from (pn| 

Lo(a^,i,C),/)(f,t,C) =0, (43) 
Lo(x, C) = t, s, C)L{x, t, C)k{x, t, s, C) = 

= -ihdt + {x{t, c), ihv)] + Sit, c) + {p{t, c),x)- {P{t, c),x{t, c))+ 
+ i((z + z{t, c)),'HUm + z{t, C))) + {ns), [z + Z{t, C))) + 

+ Z(t, C)), W,,{t)iz + Z{t, C))) + {{z + Z{t, C)),W,^{t)Z{t, C)) + 
+i(^(t, C), W^^^(t)Z(t, C)) + i Sp [W^^(t)A(2)(t, O] 



Putting 



where 



E 

5(i,C) = l[{P{t,C)j{t))-H^{t,C)}dt, (44) 



i/^(t,C) = -(Z(i,C),[H,,(t) + MW^,,(t) + 2W^,„(t) + W^^^(t))]Z(t,C)) 
+ (H,(<), Z(t, C)) + ^f,Sp[W^Ut)^^^\t, C)], 
and taking into account (jl7p . we obtain an equation for the function (j}{x,t, C): 

Loix,t)<j){x,t,C) = 0, io(f,<) = -«^5t + + }c^{z,W,,{t)z). 



(45) 



Therefore, the operator Lo{x, t) and the function (/)(x, C) do not depend on the constants C and 
equation p6p for the fundamental intertwining operator can be written as 



k{x, t, s, C')Lo{x, t)k-^{x, t, s, C')t>{t, s, C, C) = 
= Vit, s, C, C)i'(f, t, s, C)Loix, t)k-^{x, t, s, C), 

V{t,s, C,C') 



Hence, 



V{t, s, C, C) = i^-i(x, t, s, C')V{t, s, C, C)k{x, t, s, C), 
where T>{t, s, C, C) is a symmetry operator of equation pSt . i.e. 



(46) 



(47) 



Lo(f,i),X'(t,s,C',C) 



= 0, 



I?(t,s,C',C) =I?o(C',C). 



Here 



(48) 



150(0', C) = k{x,t,s,C')k-^{x,t,s,C) =exp(^<5S'(t,C,C')| 

t=s In J 

{5Xo{C, C), V) - ^{SPoiC, C'),x)] exp {^(5Xo(C, C), <5Po(C, C'))}. (49) 



X exp 



SSit,C,C') = S{t,C')^Sit,C), SZo{C,C') = (5Po(C,C'),<5Xo(C,C')), (50) 
6Xo{C, C) = Xo(C') - Xo{C), 6Po{C, C) = Po(C') - Po(C). 

The solution of the Cauchy problem (|48p for the operator 2?o(C,C) can be obtained with 
standard methods (see, e.g., [221152] ) as 

V{t, s, C, C) = VoiC, C) exp {^(5Xo(C, C), <5Po(C, C'))} exp {b{t, s, C, C)}, (51) 

where 

b{t, s, C, C) = {b{t, s, C, C), Jz), (52) 
and the 2n-component vector b — b{t, s, C, C) is a solution of the Cauchy problem for the system 

b = J[n,S) + ^W,S)]b, b =SZo{C,C'). (53) 

t—s 

Then the symmetry operator A{t) for Eq. (fT5|) (or, equivalently, (jH)) can be presented as 
(|40|) . where the intertwining operator I?(t, s,C',C) is defined by (|47p and ([51]) . and B{t,C) is a 
symmetry operator of L AE ([20]) . 

Using the explicit form ([TT]) . ([5T]) of the intertwining operator I?(i, s, C, C) and the operator 
K{x, t, s, C) from (liTI) . we have 

*A(a',i) - (i(t)*)(f,t) =exp{^[5A(i) + (PA(t),5'-X^(i))]}^(f + X(t)-X^(t),t) x 

X exp { - [5(t) + (Pit), X - X(i))] } ^'(a- + X(i) - XAit),t), (54) 

where 

i3(f,t) = B{t,C[^]{t)). 

We note that expression (|54p for the symmetry operators is not simple and requires further 
analysis, but other forms of symmetry operators for the GPE are unknown. 

To obtain more simple examples of symmetry operators in explicit form, we consider the ID 
case of equations ([T]), dHJ). 

4 Symmetry operators in ID case 

The reduced ID GPE dU) reads: 

F{^)ix,t) ^ {-ihdt + Hq + ><%{^){t)}^{x,t) = 0, (55) 

where we have used the notation 

+ 00 

- 1 ^ 1 f 

— oo 

p = —ihd/dx; a, b, and c are the real parameters of the nonlocal operator V['i>]; pi, tr, p are the 
parameters of the linear operator Hq] x, ?/ e R^. 

The Hamilton-Ehrenfest system ([TT]) for the first-order moments becomes [54] 

ip=-pp-aox, ^^^^ 
\x = fj,p + px, 

(2) (2) 

and for the second-order moments with — Aj^2 have 

AP,'.-2,Af)-2aAr), 

A^^^=pA[\^-aAg\ (58) 
Ag^ =2pAi'^ +2pA^^l 
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where 

<7o — (7 + yt(a + 5), (7 = (T + Ka. 

We introduce the notation 



\J oqii - f7 = yja^i- (59) 

and assume that fj^ = coM" > 0- This imphes that the ions do not leave the Paul trap. Indeed, 
in this case, the general solution of system (|57p is 



X{t, C) = Ci sin m + C2 cos Vlt, 

1 - - 1 - - (60) 



and all solutions of system (1571) are localized. 

Assume that the wave packets that describe the evolution of ions by Eq. (|55|) do not spread. 
This takes place if VL^ = ap — > Q. 

For system ([55]) we have 

A^2^ (i, C) = C3 sin 2nt + C4 cos 2fit + C5, 

A^i^(t, C) = -(llCg - PC4) cos 2m - -{nCi + pCz) sin 2m - -C5, 
M M M 

/o^ 1 (61) 

Aii^(t,C) = — ((p^ - ^f)Cj.+2pnCA) sin2m+ 
M 



^((p2 - ^"^^Ci - 2pl7C3) cos2m+ -( 
/i^ ^ ' p 



and all solutions of system (|58l) are also locahzed. Here, C = (Ci, . . . , C5), and C;, Z = 1, 5, are 
arbitrary integration constants. 

The associated linear equation ([20)1 in the ID case is 

L(a;,i,C)$ = 0, 

L(., C) ^ -^fta. + ^ + ^ + ^^^^ + (62) 
+xbxX{t,C) + [x2(t,C) + A^2^(t,C)]. 

We can immediately verify that for the associated linear equation (j62p we can construct the 
following set of symmetry operators linear in x and p: 

a{t, C)) = -1= [C{t){p - P(i, C)) - B{t){x - X{t, C))] , (63) 



a+it, C)) ^ -1= [C*(i)(p - P(i, C)) - B*it){x - X(t, C))] . (64) 



(65) 



Here the functions B{t) and C{t) are solutions of the linear Hamiltonian system 

(B^-pB-aC, 
[C = fiB + pC. 

The Cauchy matrix X{t) for the system ([^5)1 can easily be found as 



(cos rJi — — sin fit —(ri^+p^)smQt\ 
„ " p h ^it) =12X2. (66) 

^ sin nt cos + ^ sin fit J *=° 

The set of solutions normalized by the condition 

B{t)C* (t) - C{t)B* {t) = 2i (67) 
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can be written as 



Equation (|67p results in the following commutation relations for the symmetry operators (|63p and 

[a(i,C),a+(i,C)] =1. 
For the function (f) given by (jH]) in ID case, we obtain 



C) = K{x, t, s, C)(/)(x, i), 
K{x,t, C) = exp[-X(t, C)aj exp |^[S'(t, C) + P(t, C)x]}, (69) 

where 

t 

S{t,C) = J[p{t,C)X{t)-H^{t,C)}dt, (70) 



From (liSl) we find 



Jf^(t, C) = c) + -X2(t, C) [ao + ~>c{h + c)] + pP(i, C)X(i, C) + i^-Ag'(t, C) 



L„(.,t)0 = O, Mx,t) = -z;.ft + ^+(^±|^ + ^(^^. (71) 

Then the symmetry operator A{t) (|47)) for Eq. ((55|) can be presented as 

(i(t)vl/)(x,t) =P(t,C[a^],C[vI.](t))S(t,C[vI/](i))vI/(x,t), (72) 

where B{t, C) is the symmetry operator of the associated linear equation (j62p . The intertwining 
operator V{t, C, C) can be presented as 

V{t, C, C) = -ft:"^(a;, C')V{t, C, C)7^(a;, C), (73) 
T){t, C, C) = exp I ^ [S'(i, C) - ^(t, C)] } exp |S(t, C, C) } x 

X exp {^^1^ - C^) - p{C'^ - C^)) }, (74) 

where 

b(t, C, C) = h^{t, C C)p - 6p(i, C, C)a; = (6(i), f), (75) 
and the vector h{t) is defined by 

ftr/ c c\ - (^^^^^ - ^(^2 - c^) 

^ ' ' ^&x(i,C',C)y " i^l7(Ci-Ci)-p(C^-C72) 



The matrix A'(i) is given by (j66p . 

The symmetry operator A(t) from (|72p of nonlinear equation ([55]) has the structure of a linear 
pseudodifferential operator whose parameters are functionals of the function on which the operator 
acts. 

Therefore, the explicit form of the operator A{t) is determined not only by the symmetry 
operator of the associated linear equation B{t,C), but also by the function 'ii{x,t). Note that 
for some values of the parameters (more precisely, for the function \I'(x,t) that define them) the 
pseudodifferential operator becomes a differential one. 

We set 

B{t, C) = B^{t, C) = ^ a+{t, C) 1/ e Z+, (76) 



where the operators a+(t, C) are defined in (|63l) . 

Substituting (1751) in (17^ . we obtain the symmetry operator, which we denote by A,^{t). 

11 



With the use of a stationary solution of the Hamilton-Ehrenfest system ([57)) . ([55)) . we simphfy 
the symmetry operators and generate a countable set of explicit solutions of the ID GPE 

A stationary solution of Eqs. ([57|) . (155)) is obtained from the general solution (jSO)) . (|CT)) if we 
take integration constants as C = 0° = (Cj*, . . . , C^), = ^ = = 0, is an arbitrary 
real constant. The stationary solution is 



Substituting (1771) in ([5^ . we obtain the associated linear equation 



L(x,t,C°)$ 0, L{x,t,C") 



ihdt 



fifp' (cr + ka)x'^ p{xp + px) _ c 



(78) 



The operator K{x^ t, C) = K{x, t, C") from (|M)) is the operator of multiplication by the function 

K{x,t,C°) = exp| 



The linear operators (|55)) and (|M)) then become 



a{t, C°) = -1= - i?(t).T] , a+(t, C°) = -i= [C*(t)p ~ ; 



(79) 



they are symmetry operators for Eq. (|78| . 
The function 



M^.t,C^) = (i)''\)!)'''exp{ - A^.^ _ ^^-1exp{ - '^nt-^^kcClt} (80) 



is immediately verified to be solution of Eq. (|78)) . 

Upon direct substitution, we see that for the function ()80p . Eqs. (P^ . ([^5)) . which determine 
the functionals C[^](t), become 

X(0,C) = xv, = 0, F(0, C) = = 0, 



a('2^(0,C) = (A 



(2)^ 
22 



|c^(o)| 



2n' 



A(2i)(t,C) = (AilO^ = ^|S(0)r = 



.(2) 



2 _ h{g^ + n^) 



(81) 



Af2)(0,C) = (A^^^V = J[S(0)C*(0) + S*(0)C(0)] = 
V'(a;) = $o(a;,0,C"). 

From (HJ) and ([77]) it follows that = {h(i/2Vl). From (|33|) and ([501) we find a special solution 
*o(a;,i) of the GPE ([55]): 



*o(a;,i) = *o(.T,<,C°) 



C^ = (V/20) 
1 \l/4/f7xl/4 



The symmetry operators ([75)) . (|75)) generate from (jSO)) the solutions of the associated linear 
equation ([78)) that constitute a Fock basis in the space L2{R: 

$,(a;, t, C"') = B,(i, C°>o(x, C°') = 4=f ^'")\ ^"0 = 



W-^a; j$o(a;,i,C°')exp| -ir^z/tj, z^ e 



(83) 
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where H^{C,) are the Hermite polynomials [SS] 

i?.(C)=(2C-|)^i. 

The operator (|T7|) . intertwining the operators L(a;, t, C") and t, C'^') of the form ([751) . reads 

P(t,CO,C"') =exp{^i.c[C°-C°']i}. (84) 

Equations ([231) that determine the functionals C[^^](t) for the fimctions (|83p are written as 
X(t,C) =0, P(0,C)=p^^ =0, 

a(^2)(0,C) = (A(^2^)^<. = f (2^ + l)|C(0)|' = 1^(2^ + 1), 

Af,'(0,C) = (Af,)),„ = ^|B(0)|^(2.+ 1) = ^^^y^(2- + l), (85) 

Af2^(0, C) = (A(f )^„ = J [S(0)C*(0) + S*(0)C(0)] (2^ + 1) = -g(2^ + 1), 
^a(x) = $,(a;,0,C°') = B,(0,C°V(a;)- 

Here we have used the standard properties of Hermite polynomials [5S] . Taking into account 
dZIl), we find from dHS]) that C^' = {hii/^){v + 1/2). 

Then the symmetry operator A^{t) ([7^ transforms the solution '^(){x,t) of into solution 
5'i/(a;,t) of the nonlinear GPE ([55)1 according to the following relation: 



2?(t,C",e")^ a+(t.C") 



C0 = (?i^t/2f2)-C?' = (riAi/O)(i/+l/2) 



X 



Functions ([55)) constitute a countable set of particular solutions of Eq. ([55[) which are generated 
from 4'o(a;,i) by the nonlinear symmetry operator A^(t). 

Assume that C = = (C^ , C2, 0, 0, C5). This choice of the constants yields the following 
expression phase orbit ([60[) . ([6T]) : 



X (i , C| , Q ) = C! sin + cos 17i , 

p(t , c! ,c^)^- {ncl - pC] ) cos nt-- {nc^ + pcl ) sin (87) 

A(^,)(t,Ci) = Ci, A(^,)(t,Ci) = -^Q\ A(^,)(t,Ci) = ^Ci. 

Consider the action of the operator Ao(t) entering in (17^ . ([75)) on the functions ([5S|) . 
Let us write operator ([TT]) . intertwining the operators L(a:,t, C") and L{x,t,C'^), given by 
([75]). as 

P(t,C,C°) = exp[-X(i,C)a:,]exp|^[S'(i,C) - iifcC°t + P(t, C)x]}p(t, C, C°), 



where 



2?(i, C, C") = 2?o(C, C°) exp { + ^C2{nC\ - pCa) } exp {6(t, C, C°)}, (89) 

b{t, C, C") - 6,(t, C, C")p - 6p(t, C, CO)a; - {b{t),z), (90) 
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and a vector b{t) is defined by the expression 



6,(i,C,C0 



-Xt) 



nCi ~ pC2 



The matrix X{t) is given by 

Let us construct a nonlinear symmetry operator A(t,a)^ corresponding to the nonstationary 
phase orbit (|60l). (|6T1). The operator A{t,a) maps the nonstationary solution of equation ([55]) . 

((86)) . into another nonstationary solution of this equation, "^^{x^t). Consider the shift 

operator 

i3(t,C^) = i3(t,a,C^) =exp{aa+(t) -a*a(t)}, a G C, (91) 

where the operators a{t,C) and a+(t, C) are defined by expressions (|63)). (|64|. Operator ([9T|) is 
substituted in ((72|) for the symmetry operator B{t,C). 
Let us write the operator B{t, a, C^) entering in ([9T|) as 

a, C^) = exp{/3(i)p + 7(t)x} = exp { - ^/3(t)7(i)} exp{-fit)x} exp{/3(t)p}, (92) 



where 



[C*{t)a - C{t)a*] , j{t) = [B{t)a* - B*{t)a] 



2h' 



Thus, we have 



= exp { - y/3(0)7(0)} exp{7(0)a;}$,(a; - ihm,0, C"'), 
$^(x-in/3(0),0,C"') = 



(93) 



_ f I Y ( 1^ 

$o(a;-i?i/3(0),0,C"') = 

= <i>o(x,0,C°')exp{[(^^±^ 



a; -1^1/3(0) )$o(a;-i?i/3(0),0,C°'), 
ft / — p + o2 



a;;3(0) -i- 



Let us note that 
7(0) ^ 

m 

7(0) 



2/1 



[-B*(0)a + B(0)a*] 
[C*{0)a-C{0)a*] = 



1 



/2ri 
1 f^i 



a H : — a' 



/rip 



/Qfi 



/3'(0)]}. 
V2 



\/hfip 



[pa2 + ^ai] , 



/3(0) = ^ [ - S*(0)a + B{0)a* + _^±^(C*(0)a - C(0)a* 
P V2h L M 



^ [{p + in)a + (-P + in)a*) + {-p + in){a - a*)] 



/2hp 



y/Tinp 

Here, ai = Re a and a2 = Ima. Similarly, we have 



ih 
T 



/?(0)7(0) + ^±^/32(0) 
p 



ih2iV^a 1 r 1 r 2 I 121 • 



Substituting dHU) in (1231), we obtain 

V'a(2:) =B(0,a,Ci)$,(a:,0,Ci) = 

2iVf^ 



cxp|— ^|exp|- 



/ 1 \l/4/fl\l/4 (- 4 /-p + l\l\ -] 



zax 

/2hp J ^jv 



/2V 

IT 



-q;2 



(0 + iilN 
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(95) 



From (|95p . in particular, it follows that 



I n 

■nhii 



cxp 



{ h^i{ 



■OL2 



2-1 1 



n r 
hjj, 



■a2 



(96) 



Similar to ([55)) . we write Eqs. (^5)) determining the fmietionals C[?/'a](i) for functions ([55]) as 



1 



V2hn 



-Qfl, 



a(1)(0,C)) = (A^i^)^ = |(2^+ 1)|C(0)|' = ^{2u + l) 



Afi)(0,C)) = (AilO 



(2)^ 



|S(0)| (2z/+l) 



2r2/x 



(21^+1), 



(97) 



Ag)(0,C) = (A(f)^ = -[i?(0)C*(0)+i?*(0)C(0)](2^ + l) 
From equations (pO]) and (p7)) . in view of p2p . we obtain 



hg 

2n 



(2i. + l). 



217 



(98) 



(2zy + l). 



Then the nonlinear symmetry operator for the nonlinear GPE (1551) . A{t, a), determined by (17^ 
and ([91]) . transforms the solution '^^{x,t) ([SG)) into a nonstationary solution ^i,(x,i): 



= (i(t, a)^'^)(a;, t) = 2?(t,C\C° 
/ 1 x''/ 1 \i/4/riNi/4 



Cf=Ci = (?iA'/f2)(!^+l/2), 



^'i.(a;,t) 



2h 2h fi 



Ax' )■ X 



exp { ^ [5(t, , + P{t, Cl (a) , {a)) Ax] } , 



(99) 



which is localized around a phase orbit (a), (a)), v'C(t, Ci(a), C2 (a))). Here Ax = 

X — X{t,Cl{a),C2{ct)) with the constants Cl{a), C2(a) determined by equation (|98|) . and a 
function C|(a), Ca^a)) is determined by ^ where Ci = C|(a), C2 = C^ia), C3 = 0, 
C4 = 0, CO = (;iA'/^^)(i^ + l/2). 

In a linear case {x = 0), operators ([99)) with a G C constitute a representation of the 
Heisenberg-Weyl group [22123] . The functions '^,^{x, t, a) determined by ([99]) minimize the Schrodinger 
uncertainty relation for v — [56], and, hence, they describe squeezed coherent states [57] . 



5 Discussion 



Direct calculation of symmetry operators for nonlinear equations is, as a rule, a severe problem 
because of the nonlinearity and complexity of the determining equations [53] . 

However, for nearly linear equations [48] . a wide class of symmetry operators can be constructed 
by solving linear determining equations for operators of this type much as symmetry operators are 
found for linear PDEs. 

We have illustrated this situation by the example of the generalized multidimensional Gross- 
Pitaevskii equation ([T|). The formalism of semiclassical asymptotics leads to the semiclassically 
reduced GPE ([6]) (or ([TS]) ) which belongs to the class of nearly linear equations. Note that the 
solutions of the GPE are found in a special class of functions decreasing at infinity [47] . 

The reduced GPE is quadratic in the space coordinates and derivatives and contains a nonlocal 
term of special form. 
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In constructing the symmetry operators for the reduced Gross-Pitaevskii equation (|15p . we use 
the fact that this equation can be associated with the hnear equation (PO)) . 

The symmetry operator A{t) (^0]) of the reduced GPE (fTSjl has the structure of a hnear pseudo- 
differential operator with coefhcients C depending on the function 'if on which the operator acts. 
The operator A is given in terms of the hnear intertwining operator V and of the symmetry 
operators of the associated hnear equation ([20| . The dependence of the coefficients C on ^ arises 
from the algebraic condition ([32]) . and therefore the operator A{t) is nonlinear. This is the key 
point of the approach. 

The ID examples considered show that for a special choice of the parameters C, we can con- 
struct the symmetry operators and generate the families of solutions of the nonlinear equation (|55p 
written in explicit form. 

In conclusion, we note that the class of nearly linear equations can be considered as an approxi- 
mation of highly nonlinear equations in the context of the semiclassical approximation formalism. 

The further development of the study of symmetry operators is seen as a generalization of the 
approach for integro-differential GPE of more general form and for systems of equations of this 
type. 
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